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Finite Temperature QCD Interfaces Out of Equilibrium 
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The properties of interfaces in non-equilibrium situations are studied by constructing a density matrix with a 
space-dependent temperature. The temperature gradient gives rise to new terms in the equation for the order 
parameter. Surface terms induced in effective actions by abrupt temperature changes provide a natural theoretical 
framework for understanding the occurence of both continuous and discontinuous behavior in the order parameter. 
Monte Carlo simulation of pure QCD shows both kinds of interfacial behavior. Perturbation theory predicts a 
universal profile in the high temperature phase, which can be tested by Monte Carlo simulation. 



1. INTRODUCTION 

Lattice gauge theory is a well understood the- 
oretical approach to the equilibrium properties 
of finite temperature QCD, especially to bulk 
thermodynamic properties. Experimentally, fi- 
nite temperature effects in QCD are of great inter- 
est for studies of the early universe and heavy ion 
collisions, both non-equilibrium situations. Much 
of the theoretical uncertainty is associated with 
the interface between hadronic matter and quark- 
gluon matter. 

It is well established that computer simulations 
can give information about equilibrium interfaces, 
including the surface tension and the interfacial 
width. In principle, information about transport 
coefficients can be obtained from equilibrium cor- 
relation functions, although there are technical 
difficulties. 

An alternative approach is to examine the in- 
terfacial behavior given basic information about 
the system at a given time. For example, we 
might know that the average energy density in 
a region of space has some given value, while the 
rest of space has another. Given this information, 
it is natural to contruct a density matrix which 
maximizes the entropy subject to that constraint 
on the energy. 

The local equilibrium formulation of Zubarev 
|l) constructs the time evolution of the density 



matrix as successive changes in parameter fields 
like the local temperature and the local velocity. 
The time evolution is understood as motion from 
one local equibrium state to another. In the case 
of a system at rest, the density matrix has the 
form 



P = -^exp{- I (i{x)n{x)) 



(1) 
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where H is the energy density operator. Such 
a temperature profile can be imposed in a lattice 
simulation by noting that the inverse temperature 
(3 is given by Nta, where Nt is the temporal size 
of the lattice in units of the lattice spacing a. The 
length of a in physical units is specified by making 
the bare coupling constants space-dependent. In 
principle, this can be done in QCD in such a way 
that the spatial lattice spacing remains con- 
stant while the temporal lattice spacing at varies 
1^ ; in practice this is unnecessary in several cases 
of interest. In applications to cosmology, tem- 
perature differences may be small. This will cer- 
tainly be the case if the nucleation rate is large 
and only small undercoolings are possible. In a 
first approximation to the environment in a heavy 
ion collision, the temperature inside the region of 
interest can be taken as roughly constant at a 
high temperature, and the region outside as very 
cold. Monte Carlo studies |^ have shown that the 
Polyakov loop interfacial profile is not sensitive to 
the outer region temperature, if it is sufficiently 
cold. 
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2. BOUNDARY CONDITIONS AND 
SURFACE TERMS 

At the equilibrium point of a first-order phase 
transtion, a stable interface between two phase 
can exist. The interfacial profile for any order 
parameter is described by a one-dimensional kink 
solution which interpolates between degenerate 
minima of the effective potential. Out of equi- 
librium, the minima are no longer degenerate. In 
the traditional image of the ball rolling in an in- 
verted potential, the ball starts to roll down from 
one hilltop, and overshoots the other. Thus, it 
is necessary to understand dissipative effects to 
study the properties of a non-equilibrium inter- 
face. As shown below, a temperature gradient is 
associated with dissipative mechanisms that allow 
the existence of an interface in non- equilibrium 
situations. An abrupt change in the temperature 
will give rise to a localized dissipation at the in- 
terface. In this way, surface terms are created in 
effective actions. 

2.1. Semiclassical Argument 

Consider a simple scalar field theory, and imag- 
ine that the temperature T(z) is a given function, 
slowly varying over a mean free path of all par- 
ticles. Semiclassical arguments show that an ef- 
fective action for the one one-dimensional field 
theory is given by 



sity across the barrier to the dissipative term: 
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where (j> is the order parameter, A is the cross- 
sectional area transverse to z and Vef / is the finite 
temperature effective potential. If 14// has two 
minima, the equation for the interface is given by 
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No stable interface can form at a constant tem- 
perature other than the critical temperature, 
where the two minima are degenerate. However, 
in the presence of a temperature gradient, this 
can occur, due to the dissipative term in the ef- 
fective action. There is a kind of virial theorem, 
which relates the change in potential energy den- 
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where AVe// is the difference in the effective po- 
tential between z — +oo and z = — oo. Dissipa- 
tion takes place only in those regions where the 
temperature is changing. 

The semiclassical analysis indicates that a lo- 
calized temperature gradient will produce a local- 
ized dissipation of energy. It is easier, however, 
to analyize the effects of a sharp temperature in- 
terface from a lattice point of view. This corre- 
sponds precisely to the role of surface terms in 
interfacial problems in condensed matter physics. 

2.2. Renormalization Group 

The renormalization group provide a simple ex- 
planation for the occurence of surface terms 
Suppose a sharp planar interface at z=0 is intro- 
duced into some lattice model, with couplings J\ 
for z < and J2 for z > 0. Any real space renor- 
malization group transformation applied to this 
model will yield new bulk couplings J( and 
away from the interface. At the interface, a new 
set of couplings will be generated naturally. In 
this manner, even a microscopically sharp inter- 
face will naturally generate surface terms in an 
effective action. 

2.3. Interfacial discontinuity 

At the critical temperature of a first order 
phase transition, the order parameter changes 
continuously across an interface. Out of equi- 
librium, an expanding bubble will evolve toward 
a sharp temperature discontinuity |^. The or- 
der parameter may or may not exhibit discon- 
tinuous behavior. It is more common in con- 
densed matter studies of interfacial phenomena 
to think of bulk first order transitions as giving 
rise to first-order interfacial transitions; the pos- 
sibility of second-order interfacial behavior at a 
first-order bulk transition was first pointed out in 
1982 by Lipowsky §. 

2.4. The interfacial profile in strong and 
weak coupling 

The strong coupling limit is a natural limit in 
which to study the interfacial behavior of pure 
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Figure 1. Polyakov loop profile 



QCD, since the system simplifies to a system of 
Polyakov loop spins. Previous work using a mean 
field theory applied to an SU(3) spin model de- 
rived from strong coupling QCD shows that in a 
slab or semi-infinite geometry the surface behav- 
ior is first order. This is in reasonable agreement 
with Monte Carlo data |Q. 

Perturbative techniques can also be applied, us- 
ing a variant of the effective action developed for 
order-order interfaces in high-temperature QCD 
j?!^. This gives the interfacial profile for the 
Polyakov loop on the high-temperature side of 
the deconfinement-confinement interface. The in- 
terfacial profile is a universal function of gTz. 
However, these methods are only useful in the 
high temperature regime, and one cannot obtain 
a complete solution without an understanding 
of the low temperature confining regime. Sur- 
face terms can be added to the effective action 
to parametrize our ignorance; depending on the 
strength of these terms, the profile may go contin- 
uously to zero, or jump discontinuously to zero. 

3. MONTE CARLO RESULTS 

It is easy to simulate pure QCD at two differ- 
ent temperatures, one above and one below the 



deconfinement transition. The gas is at rest on 
both sides of the interface, and the location of 
the wall is determined by temperature profile. In 
the actual time evolution of the system, the in- 
terface is moving, in such a way that total energy 
is conserved. 

Figure 1 shows the Polyakov loop profiles for 
simulations of hot, deconfincd spherical droplets 
of radius R — 7 at Nt — 2. A spherical geometry 
was used, since it provides a detailed interfacial 
profile. The spherical surface introduces correc- 
tions to the droplet free energy which are propor- 
tional to the radius, but this disadvantage is mild 
compared to the rather coarse profile provided 
by a planar interface. The coupling constant Pout 
outside the bubble was set at 1.0; no significant 
change in the profile was seen at /3out = 3.0, al- 
though boundary values changed slightly. It was 
necessary to go to rather large values of /3i„ , cor- 
responding to high temperatures, to obtain an ap- 
parent sharp discontinuity in the Polyakov loop. 
As shown in figure 2, a mild discontinuity occurs 
when Pin = 5.6 at Nt = 2. Given that Pc — 5.69 
at A^t = 4. this correspond to a temperature of 
almost 2Tc. A fit to the continuum profile |^] 
for Pin = 5.6 gives an excellent fit; however, the 
fitted value for gT is approximately 60% higher 
than the perturbative value. A test of the per- 
turbative prediction for the profile at A^ = 4 or 
larger is indicated. 
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